DREDGING  RESEARCH  PROGRAM 

TECHNICAL  REPORT  DRP-92-6 

ADCIRC;  AN  ADVANCED  THREE-DIMENSIONAL 
CIRCULATION  MODEL  FOR  SHELVES, 
COASTS,  AND  ESTUARIES 

Report  6 

DEVELOPMENT  OF  A  TIDAL  CONSTITUENT 
DATABASE  FOR  THE  EASTERN  NORTH  PACIFIC 


3TIC 

LECTE 

lY  1  0  1995 


D 


J.  L.  Hench,  R.  A.  Luettich,  Jr. 

Institute  of  Marine  Sciences 
University  of  North  Carolina  at  Chapel  Hill 
Morehead  City,  North  Carolina  27514 

J.  J.  Westerink 


1 9950509  064 


Department  of  Civil  Engineering  and  Geological  Sciences 
University  of  Notre  Dame 
Notre  Dame,  Indiana  46556 

N.  W.  Scheffner 

Waterways  Experiment  Station,  Corps  of  Engineers 
3909  Halls  Ferry  Road,  Vicksburg,  Mississippi  39180-6199 


December  1 994 
Report  6  of  a  Series 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


Prepared  for  DEPARTMENT  OF  THE  ARMY 
U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


The  Dredging  Research  Program  (DRP)  is  a  seven-year  program  of  the  US  Army  Corps  of  Engineers. 
DRP  research  is  managed  in  these  five  technical  areas: 

Area  1  -  Analysis  of  Dredged  Material  Placed  in  Open  Waters 

Area  2  -  Material  Properties  Related  to  Navigation  and  Dredging 

Area  3  -  Dredge  Plant  Equipment  and  Systems  Processes 

Area  4  -  Vessel  Positioning,  Survey  Controls,  and  Dredge  Monitoring  Systems 

Area  5  -  Management  of  Dredging  Projects 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such 
commercial  products. 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Dredging  Research  Program 
Report  Summary 


ADCIRC:  An  Advanced  Three-Dimensional  Circulation  Model  for  Shelves,  Coasts, 
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ISSUE:  Dredged-material  disposal  sites  lo¬ 
cated  in  open  water  are  classified  as  either  dis¬ 
persive  or  nondispersive  depending  on 
whether  local  water  velocities  are  strong 
enough  to  erode  and  transport  dredged  mate¬ 
rial  from  the  deposited  mound.  The  Corps 
needs  the  capability  to  predict  stability  of  the 
mound  and  long-term  migration  patterns  of 
eroded  material  to  (1)  identify  acceptable 
disposal-site  locations,  and  (2)  provide  a  quan¬ 
titative  approach  for  gaining  site-designation 
approval. 

RESEARCH:  The  overall  work-unit  objec¬ 
tive  is  development  of  a  systematic  approach 
for  predicting  the  dispersion  characteristics  of 
a  specific  open-water  disposal  site.  This  ob¬ 
jective  includes  the  following  goals: 

•  Identify  realistic  wind-,  wave-,  tide-, 
and  storm-generated  velocity  boundary 
conditions. 

•  Develop  numerical  models  capable  of 
simulating  dispersion  characteristics  of 
dredged-material  mounds  for  periods  of 
time  in  excess  of  one  year. 

•  Provide  site-designation  technology  to 
field  engineers  as  a  tool  in  site  identifica¬ 
tion  and  designation. 


This  study  feports  the  tidal-constituent 
database  developed  for  locations  along  the 
eastern  North  Pacific  coast  of  the  United 
States,  which  is  a  part  of  the  first  two  goals 
listed  above. 

SUMMARY:  The  numerical  model 
ADCIRC-2DDI  (a  two-dimensional  depth- 
averaged  barotropic  hydrodynamic  model) 
was  applied  to  the  eastern  North  Pacific  to  de¬ 
velop  a  tidal  constituent  database.  The  report 
includes  the  specification  of  a  geometrically 
and  hydrodynamically  simple  open  boundary, 
the  use  of  large  domains,  and  the  advantages 
of  using  a  graded  finite  element  grid  to  selec¬ 
tively  resolve  flow  features  of  interest. 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Library,  telephone 
number  (601)  634-2355.  National  Technical 
Information  Service  (NTIS)  report  numbers 
may  be  requested  from  WES  Librarians. 

To  purchase  a  copy  of  the  report,  call  NTIS  at 
(703)  487-4780. 
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Summary 


This  report  describes  the  application  of  ADCIRC-2DDI,  a  two-dimensional, 
depth-averaged,  banotropic  hydrodynamic  model,  to  the  eastern  North  Pacific 
in  order  to  develop  a  tidal  constituent  database.  Issues  that  are  emphasized  in 
the  development  of  the  Eastern  North  Pacific  Tidal  (ENPACT)  model  include 
the  specification  of  a  geometrically  and  hydrodynamically  simple  open  bound¬ 
ary,  the  use  of  large  domains,  and  the  advantages  of  using  a  graded  finite 
element  grid  to  selectively  resolve  flow  features  of  interest 

The  ENPACT  model  is  coupled  to  a  finite  element  global  ocean  model  at 
the  open  boundary  and  forced  with  five  diurnal  and  semidiurnal  astronomical 
tidal  constituents  (K^  Oj,  M2,  S2,  and  N2).  Tidal  potential  forcing  is  applied 
to  the  interior  of  the  domain  for  the  same  constituents.  The  structures  of  the 
various  tides  are  examined  and  the  results  of  the  tidal  simulations  are  com¬ 
pared  with  field  data  at  72  fixed  tidal  elevation  stations. 


1  Introduction 


This  report  presents  the  initial  results  from  a  numerical  tidal  model  of  the 
eastern  North  Pacific  which  encompasses  the  coastal  ocean  as  well  as  the  deep 
open  ocean.  The  simulations  were  performed  using  ADCIRC-2DD1,  a  two- 
dimensional  depth-integrated  barotropic  finite  element  hydrodynamic  code. 

The  Eastern  North  Pacific  Tidal  (ENPACT)  model  is  intended  to  provide  a 
bridge  between  global  ocean  models  which  are  adequate  for  the  deep  ocean  but 
poorly  resolve  shelf  and  nearshore  areas,  and  small  domain  coastal  models 
(e.g.,  Flather  1987,  Forman  et  al.  1993)  which  adequately  resolve  the  nearshore 
features  but  can  suffer  from  poor  open  boundary  conditions,  particularly  from 
the  higher  harmonics  which  are  generated  on  the  shelf.  The  results  of  this 
model  are  intended  to  provide  a  tidal  constituent  database  from  which  accurate 
boimdary  conditions  can  be  extracted  to  drive  smaller  localized  models  with 
higher  resolution.  In  earlier  work,  a  similar  database  was  generated  for  the 
entire  western  North  Atlantic  using  the  same  code  (see  Westerink,  Luettich, 
and  Scheffher  1993). 

Issues  emphasized  in  this  report  include  the  use  of  a  geometrically  and 
hydrodynamicaUy  simple  open  boundary,  the  use  of  a  large  computational 
domain  in  order  to  accurately  propagate  tides  from  the  deep  ocean  onto  the 
shelf,  and  the  use  of  a  graded  finite  element  grid  to  reduce  computational  over¬ 
head.  Model  results  are  compared  with  observational  data  at  72  fixed- 
elevation  recording  stations  distributed  throughout  the  domain. 

The  size  of  the  ENPACT  model  domain  is  nearly  2  x  10^  sq  km  and  would 
be  computationally  prohibitive  to  discretize  in  a  uniform  manner  at  the  resolu¬ 
tion  required  to  resolve  the  rapidly  varying  flow  features  on  the  shelf.  Much 
of  the  domain  lies  in  the  deep  ocean,  which  requires  only  relatively  coarse 
resolution  since  the  tides  vary  so  slowly  there.  Therefore,  a  graded  finite  mesh 
was  designed  to  minimize  the  size  of  the  discrete  problem,  while  stUl  resolving 
the  flow  features  of  interest. 

The  remainder  of  this  report  discusses  the  application  of  ADCIRC-2DDI  to 
simulate  astronomical  tides  in  the  eastern  North  Pacific.  Chapter  2  gives  a 
brief  summary  of  the  formulation  and  numerical  discretization  of  the  model. 
Chapter  3  describes  the  computational  domain  and  the  modehng  strategy. 
Chapter  4  details  the  model  runs.  In  Chapter  5,  comparisons  between  modeled 


Chapter  1  Introduction 
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and  observed  tidal  elevations  are  made.  Chapters  5  and  6  provide  discussion 
and  conclusions  for  this  study.  Finally,  Appendix  A  outlines  procedures  for 
resynthesizing  harmonic  constituents  for  model  results  and  field  data. 
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Chapter  1  Introduction 


2  Governing  Equations  and 
Numerical  Discretization 


The  computations  described  in  this  report  were  performed  using 
ADCIRC-2DDI,  the  depth-integrated  version  of  a  group  of  two-  and  three- 
dimensional  hydrodynamic  codes  named  ADCIRC.  A  full  treatment  of  the 
governing  equations  and  the  numerical  techniques  used  in  ADCIRC  is  given  in 
Luettich,  Westerink,  and  Schefftier  (1992).  The  following  briefly  outlines  their 
derivation  and  numerical  discretization.  ADCIRC-2DDI  uses  the  depth- 
integrated  equations  of  mass  and  momentum  conservation,  subject  to  the 
incompressibility,  Boussinesq,  and  hydrostatic  pressure  approximations.  Using 
the  standard  quadratic  parameterization  for  bottom  stress  and  neglecting  baro- 
clinic  terms  and  lateral  diffusion/dispersion  effects  leads  to  the  following  set  of 
conservation  statements  in  primitive  non-conservative  form  expressed  in  a 
spherical  coordinate  system  (Flather  1988,  Kolar  et  al.  1994): 


ac  ^  1  [au/f  ^  acw  cos  (p)]  ^ 

"a?  R  cos  9[  ax,  atp  J 

+  -ftancPn  ^f]v  = 

a?  R  cos  9  ax  R  d(p  y  R  J 


R  cos  9  ^  Po 


_L  +  g(C  -  n)  h -11  -  T.u 


(1) 


(2) 


dt  R  COS  9  ax  R  d<f> 


R 


U  = 


1  ^  fi 

R 


g(C  -  n) 


.  Jfl  -  T.V 


(3) 
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where 

C  =  free  surface  elevation  relative  to  the  geoid 
t  =  time 

R  =  radius  of  the  earth 

X,  9  =  degrees  longitude  (east  of  Greenwich  positive)  and  degrees 
latitude  (north  of  the  equator  positive) 

U,V=  depth-averaged  horizontal  velocities 

H  +  h  =  total  water  column 

h  =  bathymetric  depth  relative  to  the  geoid 

f=2Clsin<p=  Coriolis  parameter 

=  angular  speed  of  the  earth 

=  atmospheric  pressure  at  the  free  surface 

Po  =  reference  density  of  water 

g  =  acceleration  due  to  gravity 

r|  =  effective  Newtonian  equilibrium  tide  potential 

t:.<p  =  applied  free  surface  stress 


X*  = 


,  (£/2  +  y2)iy2 

- H - 


Cf  =  bottom  friction  coefficient 


A  practical  expression  for  the  effective  Newtonian  equilibrium  tide  potential 
is  given  by  Reid  (1990)  as: 
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where 


X,(p  =  degrees  longitude  and  latitude,  respectively 

Ujn  =  effective  earth  elasticity  factor  for  tidal  constituent  n  of  species  J 

Cj„  =  constant  characterizing  the  amplitude  of  tidal  constituent  n  of  species 

J 

fj„  =  time-dependent  nodal  factor 

y  =  0,1,2  =  tidal  species  (/  =  0,  declinational; ;  =  1,  diurnal; 
j  =  2,  semidiurnal) 

tQ  =  reference  time 

Tj„  =  period  of  constituent  n  of  species  J 
Vj„  =  time-dependent  astronomical  argument 
Lq  =  3sin^  (p  -  1 
Li  =  sin(2(p) 

Lj  =  cos^(((>) 

Values  for  are  presented  by  Reid  (1990).  The  value  for  the  effective 
earth  elasticity  factor,  OLj„,  is  typically  taken  as  0.69  for  all  tidal  constituents 
(Schwiderski  1980,  Hendershott  1981)  although  its  value  has  been  shown  to  be 
slightly  constituent  dependent  (Wahr  1981,  Woodworth  1990). 

To  facilitate  a  finite  element  (FE)  solution  to  Equations  1-3,  these  equations  • 
are  mapped  from  spherical  form  onto  a  rectilinear  coordinate  system  using  a 
Carte  Parallelogrammatique  (CP)  projection  (Pearson  1990): 


=  RiX  -  X,o)cos(po 

(5) 

=  R<p 

(6) 

where  Xq,  (po  =  center  point  of  the  projection.  Applying  the  CP  projection  to 
Equations  1-3  gives  the  shallow- water  equations  in  primitive  nonconservative 
form  expressed  in  the  CP  coordinate  system: 

aC  cos(po  d(UH)  ^  1  d(yH  cos  (p)  ^  ^ 

cos(p  3jc'  cos(p 
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cos(po  a  p 


-  li .  g(C  -  n) 

COS(p  ax'  Po  Po^ 


av  cos  (pn  av  dV  f  tanm  ^ 

^  + _ ZlUZL  +  -  '^'rn  +f  u 

W  cos  9  ax'  dj/  1  /?  J 

r  1 

=  -^  £i.s(C-n) 

9}''  [Po  J  Po^ 

Utilizing  the  FE  method  to  resolve  the  spatial  dependence  in  the  shallow- 
water  equations  in  their  primitive  form  chronically  gives  inaccurate  solutions 
with  severe  artificial  near  2  •  Ax  modes  (Gray  1982).  However,  reformulating 
the  primitive  equations  into  a  generalized  wave  continuity  wave  equation 
(GWCE)  gives  highly  accurate,  noise-free,  FE-based  solutions  to  the  shallow- 
water  equations  (Lynch  and  Gray  1979,  Kinnmaric  1984).  The  GWCE  is 
derived  by  taking  die  time  derivative  of  the  primitive  continuity  equation  and 
spatial  derivatives  of  the  primitive  horizontal  momentum  equations.  The  spa¬ 
tially  differentiated  primitive  horizontal  momentum  equations  are  then  substi¬ 
tuted  into  the  time-differentiated  primitive  continuity  equation,  which  in  turn  is 
added  to  the  original  primitive  continuity  equation  multiplied  by  a  constant  in 
space  and  time,  Tq  (Lynch  and  Gray  1979;  Kinnmaric  1984;  Luettich, 
Westerink,  and  Scheffher  1992).  After  applying  the  CP  coordinate  projection, 
the  GWCE  is: 


^  a  8C 

°  dr  COS9  dx'  9r 


U  -  -  vhUL 

COS9  dx'  dy' 


i£L2t7  ./Vff  -  _ilfi  *  s(;  -  ti) 

R  j  COS9  a;c'|_Po 
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Po  ay'  C0S9  dx' 
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The  GWCE  (Equation  10)  is  solved  in  conjunction  with  the  primitive  momen¬ 
tum  equations  in  nonconservative  form  (Equations  8  and  9). 

The  high  accuracy  of  GWCE-based  FE  solutions  is  a  result  of  their  excel¬ 
lent  numerical  amplitude  and  phase  propagation  characteristics.  In  fact  Fourier 
analysis  indicates  that  in  constant  depth  water  and  using  linear  interpolation,  a 
linear  tidal  wave  resolved  with  25  nodes  per  wavelength  is  more  than  ade¬ 
quately  resolved  over  the  range  of  Courant  numbers:  C  =  AtlAx<L0 
(Luettich,  Westerink,  and  Scheffher  1992).  Furthermore,  the  monotonic  disper¬ 
sion  behavior  of  GWCE-based  FE  solutions  avoids  generating  artificial  near 
2  •  Ax  modes,  which  plague  primitive-based  FE  solutions  (Platzman  1981, 
Foreman  1983).  The  monotonic  dispersion  behavior  of  GWCE-based  FE  solu¬ 
tions  is  very  similar  to  that  associated  with  staggered  finite  difference  solutions 
to  the  primitive  shallow- water  equations  (Westerink  and  Gray  1991).  GWCE- 
based  FE  solutions  to  the  shallow-water  equations  permit  extremely  flexible 
spatial  discretizations,  which  can  substantially  reduce  the  discrete  size  of  any 
problem  (Le  Provost  and  Vincent  1986,  Foreman  1988,  Vincent  and 
Le  Provost  1988,  Westerink  et  al.  1992). 

The  details  of  ADCIRC  and  the  implementation  of  the  GWCE-based  solu¬ 
tion  to  the  shallow-water  equations  are  described  by  Luettich,  Westerink,  and 
Scheffher  (1992).  As  with  most  GWCE-based  FE  codes,  ADCIRC  applies 
three-noded  linear  triangles  for  surface  elevation,  velocity,  and  depth.  Further¬ 
more,  the  decoupling  of  the  time-independent  and/or  tridiagonal  system 
matrices,  elimination  of  spatial  integration  procedures  during  time-stepping, 
and  full  vectorization  of  aU  major  loops  result  in  a  highly  efficient  code. 
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3  Description  of  the 
Computational  Domain 


The  ENPACT  model  domain,  shown  in  Figure  1,  includes  most  of  the 
eastern  North  Pacific,  the  Gulf  of  California,  and  the  Guatemala  and  Panama 
Basins.  The  curved  open  boundary  stretches  from  Seal  Cape  on  Unimak 
Island,  Alaska,  in  the  north,  to  Punta  Parada,  Peru,  in  the  south.  The  land 
boundary  includes  the  Pacific  coastlines  of  Alaska,  Canada,  mainland  U.  S., 
Mexico,  Guatemala,  El  Salvador,  Nicaragua,  Costa  Rica,  Panama,  Columbia, 
and  Northern  Peru. 

Although  the  focus  of  this  model  is  on  the  U.S.  west  coast  and  Alaska,  the 
larger  scale  physics  of  the  eastern  Pacific  must  also  be  considered.  A  number 
of  diurnal  and  semidiurnal  amphidromes  exist  in  the  deep  ocean  which  have  a 
substantial  effect  on  the  tides  in  the  coastal  areas.  Cartwright,  Ray,  and 
Sanchez  (1991)  established  the  approximate  locations  of  these  amphidromes 
using  Geosat  altimetry.  Figure  1  shows  the  approximate  locations  of  these 
amphidromes.  Due  west  of  southern  California  are  Mj,  S2,  and  Nj  amphi¬ 
dromes;  west  of  the  Galapagos  Islands  amphidromes  exist  for  the  Kj,  Oj,  M2, 
S2,  and  N2  constituents.  Since  spatial  variations  in  phase  are  rapid  near  amphi¬ 
dromes,  it  would  be  extremely  difficult  to  specify  accurate  open  boundary  • 
conditions  near  any  of  these  features.  For  this  reason,  the  ENPACT  model’s 
open  boundary  is  located  weU  away  from  all  of  these  amphidromes.  In  addi¬ 
tion,  the  vast  southern  extent  of  the  domain  allows  inclusion  of  large-scale 
effects  on  the  U.S.  west  coast  tides  arising  from  the  complex  tides  in  the 
Panama  Basin  and  the  Gulf  of  California.  However,  the  open  boundary  is  not 
so  far  offshore  that  the  presence  of  additional  observed  M2,  S2,  and  N2 
amphidromes  to  the  south,  and  the  Hawaiian  Islands  to  the  west,  would 
introduce  additional  complications. 

Further,  the  model  open  boundary  is  almost  entirely  in  the  deep  ocean,  well 
away  from  the  continental  shelf  and  slope  where  nonlinear  constituents  are 
important  but  difficult  to  determine.  (Where  the  open  ocean  boundary  meets 
the  land  in  the  Aleutian  Peninsula  and  in  Peru,  the  continental  shelf  and  slope 
are  very  narrow  so  the  effects  of  nonlinear  tides  in  these  areas  should  be  negli¬ 
gible.)  This  positioning  simplifies  the  specification  of  open  boundary  condi¬ 
tions  considerably.  It  also  facilitates  the  coupling  of  the  ENPACT  model  with 
global  tidal  models  which  are  most  accurate  in  the  deep  ocean.  Finally,  the 
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geometrically  simple  open  boundary  follows  a  single  smooth  continuous  arc, 
avoiding  the  complications  which  arise  when  multiple  boundaries  and  bound¬ 
aries  with  sharp  comers  are  used. 

Land  boundaries  in  the  ENPACT  model  were  derived  from  the  World  Data 
Bank  II  geophysical  database  (NGDC  1993).  This  database  provides  an  accu¬ 
rate  digital  representation  of  all  major  coastline  features  and  islands  and  is 
particularly  well-suited  for  large  scale  models. 

Model  bathymetry  is  shown  in  Figure  2.  Depths  range  from  over  6,000  m 
in  the  Aleutian  Trench  to  less  than  1  m  near  shore.  The  continental  shelf  is 
very  narrow  in  most  areas,  and  extremely  sharp  bathymetric  gradients  exist 
along  the  continental  shelf  and  slope.  AU  bathymetric  data  for  the  model  were 
obtained  from  the  ETOPO-5  database  (NGDC  1988),  which  has  a  5-min  by  5- 
min  resolution.  Depths  at  each  grid  node  were  linearly  interpolated  from  this 
database. 

By  taking  advantage  of  the  grid  flexibility  permitted  by  the  finite  element 
method,  grid  points  in  the  model  are  concentrated  in  areas  with  rapidly  varying 
bathymetry  and/or  complex  shoreline  geometry.  Conversely,  fewer  grid  points 
are  used  in  the  open  ocean  and  the  deep  Guatemala  Basin  (Figure  3).  The 
triangular  finite  element  grid  was  generated  using  ACE/gredit,  a  flexible,  inter¬ 
active,  graphical  software  package  (Turner  and  Baptista  1991).  In  areas  of  the 
grid  where  depths  exceed  3,000  m,  resolution  is  about  0.5  deg  (about  44  km). 
Between  the  1,000-m  and  3,000-m  isobaths,  grid  resolution  is  about  15  min 
(about  22  km).  North  of  Cabo  Corrientes,  Mexico,  and  above  the  1,000-m 
isobath,  grid  resolution  is  7.5  min  (about  1 1  km).  In  some  nearshore  areas 
such  as  Cook  Inlet,  Alaska,  and  Puget  Sound,  Washington,  the  grid  resolution 
is  even  finer.  Transitions  from  larger  elements  to  smaller  elements  are  smooth, 
and  all  elements  in  the  grid  passed  a  stringent  screening  for  skewness. 

Extensive  grid  convergence  studies  conducted  with  a  large-scale  model  of 
the  western  North  Atlantic  (Westerink,  Luettich,  and  Muccino  1994;  Luettich 
and  Westerink,  in  press)  suggest  that  the  ENPACT  grid  should  be  adequate  to 
resolve  most  flow  features  of  interest.  In  addition,  since  grid  resolution  in 
nearshore  areas  is  roughly  at  the  same  resolution  as  the  ETOPO-5  database, 
further  grid  refinement  in  these  areas  would  be  of  limited  benefit  without  addi¬ 
tional  bathymetric  data. 

The  final  grid  used  for  the  tidal  simulations  has  27,494  nodes.  This  is 
about  one-fourth  of  the  approximately  100,000  nodes  which  would  be  required 
if  the  domain  were  discretized  at  a  uniform  resolution  of  7.5  min. 
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4  Tidal  Simulations  with  the 
ENPACT  Model 


The  tidal  simulations  were  performed  using  the  hydrodynamic  code 
ADCIRC-2DDI  (version  25.05).  The  runs  included  the  finite  amplitude  and 
convective  acceleration  terms,  as  well  as  nonlinear  bottom  friction.  In  addi¬ 
tion,  spatially  varying  Coriolis  and  tidal  potential  terms  were  included,  both  of 
which  are  essential  when  modeling  a  domain  of  this  size.  The  bottom  friction 
coefficient  was  set  equal  to  0.003  (Chow  1959)  for  the  entire  domain.  Since 
the  runs  were  intended  to  be  entirely  predictive,  no  attempt  was  made  to  tune 
this  parameter. 

Details  of  the  model  input  requirements  and  the  use  of  its  preprocessor, 
ADCSETUP,  are  given  in  Westerink  et  al.  (1994).  Key  model  parameters  for 
the  ENPACT  model  runs  are  summarized  briefly  next.  Using  a  time-step  of 
37.5  sec,  the  minimum  and  maximum  wave  celerity-based  Courant  numbers 
(C  =  VgA  Ar/Ax)  in  the  grid  were  0.01  and  0.99,  respectively.  To  avoid  drying 
of  computational  elements,  a  minimum  depth  of  6  m  was  set  for  the  entire 
domain.  Time  weighting  factors  of  0.35,  0.30,  and  0.35  were  used  for  the 
future,  currem,  and  past  time  levels  of  the  GWCE.  The  weighting  coefficient 
for  the  primitive  portion  of  the  GWCE,  Tq,  was  set  to  0.001.  Center  point 
longimde  and  latitude  for  the  CP  projection  were:  7^  =  125°  W,  and  <P(,  = 

21°  N. 

The  open  boundary  was  forced  with  two  diurnal  astronomical  constituents, 
Ki  and  Oi;  and  with  three  semidiurnal  astronomical  constituents,  M2,  Sj,  and 
Nj.  These  boundary  conditions  were  obtained  by  inteipolating  results  from  a 
finite  element  global  ocean  model  (Le  Provost  et  al.  1994)  onto  each  of  the 
model’s  open  boundary  nodes.  In  the  deep  ocean  these  boundary  conditions 
show  good  agreement  with  the  observed  values  of  Cartwright,  Ray,  and 
Sanchez  (1991)  and  the  International  Hydrographic  Organization  (IHO)  (1990). 
In  addition,  they  j^pear  to  be  an  improvement  over  the  widely  used 
Schwiderski  (1980)  global  model,  which  was  shown  to  overpredict  many  of 
the  constituents  in  a  detailed  study  of  the  western  North  Atlantic  (Westerink, 
Luettich,  and  Muccino  1994).  The  range  of  values  used  to  force  the  open 
boundary  is  presented  in  Table  2.  Tidal  potential  forcing  was  applied  using 
the  same  constituents  as  were  used  at  the  open  boundary.  The  tidal  potential 
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constants  and  the  constituent-dependent  earth  elasticity  factors  are  given  in 
Table  3. 

The  model  was  spun  up  from  rest  using  a  smooth  hyperbolic  tangent  ramp, 
acting  on  both  the  open  boundary  and  direct  forcings  for  15  days.  The  use  of 
this  ramp  sharply  reduces  the  generation  of  short-period  gravity  and  vortex 
modes  due  to  startup  transients.  The  model  simulation  was  a  total  of  62  days 
long,  the  first  30  days  of  which  were  discarded  to  ensure  that  any  lingering 
startup  transients  had  dissipated.  A  least-squares  harmonic  analysis  was  per¬ 
formed  on  the  remaining  32  days. 

To  confirm  that  the  model  had  reached  a  state  of  dynamic  equilibrium,  the 
model  was  run  for  an  additional  33  days  (to  day  95)  and  a  separate,  indepen¬ 
dent  harmonic  analysis  conducted.  Figure  4  shows  the  amplitude  and  phase 
differences  between  the  harmonic  analysis  done  for  days  30  through  62  and  the 
harmonic  analysis  for  days  63  through  95,  for  the  M2  and  Kj  constituents. 
Amplitude  differences  for  the  Mj  tide  are  typically  less  than  0.1  cm,  and 
nowhere  in  the  domain  are  they  greater  than  0.5  cm.  Phase  differences  for  the 
Mj  tide  are  typically  less  than  1  deg  for  the  entire  domain,  and  nowhere  are 
the  differences  greater  than  2  deg  -  even  near  amphidromes.  The  Kj  tide 
shows  similar  behavior,  although  the  magnitudes  of  the  differences  are  slightly 
higher.  Kj  amplitude  differences  are  typically  1  cm,  and  phase  differences  are 
typically  less  than  1  deg.  In  Puget  Sound,  WA,  which  is  poorly  resolved  in 
this  grid,  K,  amplitude  differences  approach  5  cm  and  phase  differences 
approach  9  deg.  However,  for  the  vast  majority  of  the  domain,  there  are  no 
meaningful  differences  between  the  two  analyses.  Therefore,  the  results  of  this 
comparison  indicate  that  the  model  had  indeed  reached  dynamic  steady-state 
by  the  end  of  the  initial  32  days.  It  also  demonstrates  the  ability  of  the 
ADCIRC  code  to  run  for  over  3  months  without  generation  of  spurious  modes. 

The  computational  efficiency  of  the  model  was  very  good.  A  number  of 
techniques  have  been  implemented  in  ADCIRC  to  improve  model  performance 
including  full  vectorization  of  all  major  loops  and  the  use  of  time-independent 
and/or  tridiagonal  system  matrices.  In  addition,  ADCIRC  has  the  capability  to 
use  a  variety  of  iterative  compact  storage  solvers  from  the  ITPACKV  2D  pack¬ 
age  (Kincaid,  Oppe,  and  Young  1989).  These  iterative  solvers  have  excellent 
stability  and  accuracy  characteristics  and  are  much  more  memory-efficient  than 
conventional  banded  matrix  solvers,  especially  for  large  problems.  For  the 
27,494-node  ENPACT  grid,  the  total  CPU  time  on  a  Cray  Y-MP  8/8128  for 
the  62-day  simulation  including  a  20-frequency  harmonic  analysis  was 
19.79  hr.  On  the  same  machine,  the  maximum  memory  required  for  the  sim¬ 
ulation  was  5.87  MWords. 
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5  Comparison  Between 
Computed  and  Observed 
Tidal  Elevations  in  the 
ENPACT  Model 


The  coupled  global/ENPACT  model  results  allow  for  realistic  comparisons 
with  observed  tidal  elevation  data.  A  relatively  large  number  of  these  data  are 
available  for  the  coastal  areas  of  the  eastern  Pacific  in  harmonically  decom¬ 
posed  form  from  the  International  Hydrographic  Organization  Tidal  Constituent 
Bank  (1990).  Only  observational  data  with  record  durations  longer  than 
189  days  were  used  to  ensure  reliability  and  consistency.  However,  for  some 
places  in  the  domain  long  tidal  records  were  unavailable.  In  these  areas  (typi¬ 
cally  at  offshore  locations),  short-term  observational  data  were  used  to  improve 
the  spatial  coverage  of  the  comparisons.  Table  1  presents  the  station  locations 
and  record  durations  for  each  observational  station.  Figure  5  shows  the  posi¬ 
tion  of  each  station  within  the  model  domain. 

Model  physics  include  the  full  set  of  nonlinear  terms  in  the  shaUow-water 
equations  and  therefore  allow  constituents  to  interact  with  themselves  and  each 
ottier  to  generate  overtides  and  compound  tides.  Results  from  the  simulations 
were  harmonically  analyzed  using  half-hour  sampling  from  the  final  32  days  of 
the  simulation.  Employing  a  least-squares  method,  results  were  decomposed 
into  20  frequencies.  Although  only  the  five  astronomical  forcing  frequencies 
with  observed  data  were  compared,  fourteen  additional  overtide  and  compound 
frequencies  and  a  steady  component  were  included  in  the  harmonic  analysis  to 
more  accurately  separate  individual  constituents  from  the  tidal  signal.  The 
computations  were  intended  to  be  entirely  predictive,  so  no  calibration  or  tun¬ 
ing  procedures  were  performed  to  improve  agreement  between  model  results 
and  observed  values. 

Cotidal  charts  for  the  two  diurnal  constituents  and  the  three  semidiurnal 
constituents  are  presented  in  Figures  6  and  7.  The  diurnal  constituents  show  a 
strong  similarity  to  each  other.  Diurnal  amplitudes  show  a  steady  increase 
with  latitude  up  into  the  Gulf  of  Alaska.  The  Kj  and  Oj  tides  each  have 
amphidromes  west  of  the  Galapagos  Islands.  Much  of  die  U.S.  coast  is  nearly 
in  phase  for  the  diurnal  tides.  The  semidiumals  also  show  a  strong  similarity 
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to  each  other.  Semidiurnal  amplitudes  markedly  increase  as  these  tides  propa¬ 
gate  into  the  Gulf  of  Alaska  and  into  the  Guatemala  and  Panama  basins. 
Extremely  sharp  increases  in  amplitude  occur  in  Cook  Inlet,  AK,  and  along  the 
Canadian  coast  where  the  continental  shelf  is  relatively  wide.  The  Mj,  S2,  and 
Nj  phases  each  exhibit  twin  amphidromes  in  the  deep  ocean.  Large  stretches 
of  the  U.S.  west  coast  are  within  60  deg  of  being  in  phase  for  these 
constituents. 

The  accuracy  of  the  tidal  simulations  was  quantified  by  comparing  modeled 
values  at  the  72  stations  with  observed  values  on  a  constituent-by-constituent 
basis.  Figures  6  and  7  shows  relative  amplitude  errors  and  absolute  phase 
errors  at  each  of  the  72  stations.  To  more  readily  identify  trends  in  the  com¬ 
parisons,  the  stations  were  grouped  into  five  separate  regions.  The  station 
groupings  are  as  follows:  South  and  Central  America  (stations  1-19), 

U.S.  west  coast  (stations  20-41),  Canada  and  Alaska  (42-64),  and  Offshore 
(stations  65-72).  For  each  of  the  defined  regions  as  well  as  for  the  entire 
domain,  the  amplitude  error  for  constituent  j  for  each  of  the  five  astronomical 
constituents  was  computed  as  a  proportional  standard  deviation: 


J-,  calc-obs  _  ^ 
^j-amp  ^ 


/=1 


’ir(wi)F 


1/2 


(11) 


where 


L  =  total  number  of  elevation  stations  within  a  given  region 

(Jt/,y/)  =  location  of  an  elevation  station 

11“^  (Xi,yi)  =  computed  elevation  amplitude  for  constituent  j  at  station 
coordinates  (x,,y,) 

=  observed  elevation  amplitude  for  constituent  j  at  station 
coordinates  (xi,y^ 

The  phase  error  for  each  constituent  J  was  computed  as  an  absolute  average 
error  defined  over  a  region  as: 

-  |E  l-pfW  -  ffW,)!  ' 
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where 


(Xi,yi)  =  computed  elevation  phase  for  constituent  j  at  station 
coordinates  (x,,y,) 

^‘(xi,yi)  =  observed  elevation  phase  for  constituent  j  at  station 
coordinates  (Xi,yi) 

Amplitude  and  phase  errors  are  presented  in  Tables  4  and  5  for  the  entire 
domain  as  well  as  on  a  regional  basis  for  each  of  the  five  astronomical  constit¬ 
uents.  For  all  72  stations  in  the  domain  all  constituents  are  predicted  with  an 
average  amplitude  error  between  15.3  and  65.5  percent  and  an  average  phase 
error  between  5.0  and  66.2  deg.  In  general,  the  amplitude  errors  are  smaller  at 
the  offshore  and  South  and  Central  America  stations  where  the  tides  are  rela¬ 
tively  small.  Along  the  U.S.  west  coast  amplitude  errors  are  relatively  large, 
particularly  for  the  Sj  constituent.  Phase  errors  show  no  predominant  trends  in 
spatial  variability. 

The  distribution  of  amplitude  ratio  errors  for  each  constituent  is  presented 
in  Table  6.  In  general  the  predicted  values  for  the  diurnal  constituents  are 
closely  grouped  around  the  observed  values,  with  a  slight  bias  toward  oveipre- 
diction.  The  semidiurnal  constituents  show  quite  dissimilar  behavior  among 
themselves.  The  Mj  constituent  shows  a  somewhat  uniform  distribution 
around  the  observed  values  with  a  small  peak  at  the  mean.  The  S2  constituent 
is  in  general  grossly  over-predicted,  but  with  a  large  number  of  under¬ 
predicted  outliers  as  well.  Finally,  the  calculated  values  for  the  Nj  constituent 
tightly  grouped  around  the  observed  values  with  relatively  few  outliers. 

Table  7  gives  the  distribution  of  phase  error  for  each  constituent  Phase 
predictions  for  the  two  diurnal  constituents,  Kj  and  O,,  are  quite  good  with 
most  stations  within  ±10  deg  of  the  observed  values.  The  Mj  and  particularly 
the  Sj  tide  tend  to  lead  the  observations  by  a  large  margin.  Finally,  the  calcu¬ 
lated  Nj  shows  a  wide  distribution  of  phase  error  with  the  highest  percent  of 
stations  being  ahead  of  the  observed  phase. 

Resynthesized  surface  elevation  time  histories  derived  from  both  the  com¬ 
puted  and  observed  constituent  data  are  shown  for  representative  stations  in 
Figure  8.  Computed  time  histories  were  resynthesized  from  harmonically 
decomposed  ENPACT  model  results  for  the  five  astronomical  forcing  constitu¬ 
ents.  Observed  time  histories  were  resynthesized  in  the  same  manner.  Appen¬ 
dix  A  provides  procedures  for  resynthesis. 

These  time  histories  show  the  types  of  tides  which  occur  in  the  eastern 
North  Pacific  and  are  useful  for  making  general  comparisons  between  model 
results  and  the  observed  tides.  Along  South  and  Central  America,  the  tides  are 
mixed  semidiurnal.  The  tides  become  decreasingly  mixed  with  increasing 
latitude,  and  north  of  Washington  state  they  are  almost  completely  semidiurnal. 
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6  Discussion 


The  ability  of  a  numerical  model  to  successfully  simulate  enviroranental 
scale  flows  depends  on  the  physics  included  in  the  model,  the  numerical  accu¬ 
racy  (which  depends  on  the  algorithm  as  well  as  the  grid),  the  representation 
of  the  bathymetry  and  the  land  geometry,  and  the  precision  of  the  forcing 
functions.  Two-dimensional  barotropic  models  generally  give  comparisons 
between  model  predictions  and  measurements  within  10-percent  accuracy  for 
amplitude  and  10-percent  accuracy  for  phase.  Therefore,  it  is  assumed  that  the 
physics  incorporated  into  ADCIRC-2DDI  should  be  generally  satisfactory  for 
large-scale  tidal  simulations.  However,  it  should  be  noted  that  in  localized 
coastal  areas  where  the  water  column  is  strongly  stratified  or  where  freshwater 
inflows  are  large,  the  model  physics  may  be  not  be  adequate. 

The  high  numerical  accuracy  of  GWCE-based  finite  element  models, 
including  ADaRC-2DDI,  is  well  documented  (e.g..  Lynch  and  Gray  1979;  Le 
Provost  and  Vincent  1986;  Luettich,  Westerink,  and  Scheffher  1992).  Grid 
properties  of  the  model  arc  quite  good:  all  elements  are  close  to  equilateral  in 
shape,  adjacent  elements  are  of  similar  sizes,  and  nodal  support  is  even 
throughout  the  domaia  Earlier  grid  convergence  studies  (Westerink,  Luettich, 
and  Muccino  1994;  Luettich  and  Westerink,  in  press)  indicate  that  the 
ENPACT  model  grid  should  be  reasonably  converged  for  most  areas  within  the 
domain  for  the  K^,  Oj,  M2,  S2,  and  N2  tides. 

Inaccurate  representation  of  the  bottom  topography  in  the  model  is  likely  to 
be  a  major  source  of  error  in  the  results,  particularly  near  the  coast.  In  deep 
water.  Smith  (1993)  found  that  the  ETOPO-5  database  could  have  a  mean 
error  of  26  m.  In  nearshore  areas,  where  bathymetric  gradients  are  the  largest, 
ETOPO-5  bathymetric  data  are  the  least  accurate.  Further,  although  the 
ETOPO-5  database  has  a  uniform  resolution  of  5  min  by  5  min,  the  original 
data  finm  which  it  was  derived  have  a  highly  variable  spatial  density.  Above 
the  200-m  isobath,  ETOPO-5  soundings  are  the  most  sparse  and  least  reli¬ 
able.^  One  improvement  that  could  be  made  to  the  present  study  would  be  to 
use  more  accurate  representations  of  the  bottom  with  higher  resolution  digital 
bathymetric  data  as  they  become  available. 


*  Personal  Communication,  June  1994,  P.W.  Sloss,  NOAA-NGDC,  Boulder,  CO. 
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The  open  boundary  conditions  which  were  interpolated  from  a  global  ocean 
model  (Le  Provost  et  al.,  in  press)  are  spatially  smooth  and  in  general  show 
good  agreement  with  the  observational  data  of  Cartwright,  Ray,  and  Sanchez 
(1991)  and  the  IHO  (1990).  However,  the  coupled  global/ENPACT  model 
results  show  large  differences  from  observed  values  in  the  deep  ocean.  Off¬ 
shore  the  forced  constituent  tides  are  quite  complex  due  to  the  presence  of 
multiple  amphidromes.  Horizontal  translations  in  the  position  of  any  of  these 
amphidromes  can  have  a  basin  scale  effect  on  the  tides  (Luettich  and 
Westerink,  in  press),  although  the  areas  closest  to  the  amphidromes  will 
usually  show  Ae  largest  relative  change. 

Comparisons  with  observed  results  suggest  that  several  localized  areas  of 
the  grid  require  additional  resolution.  For  example,  on  the  Washington  shelf 
and  at  the  mouth  of  the  Strait  of  Juan  de  Fuca  model,  results  are  reasonably 
accurate.  However,  in  the  immediately  adjacent  waters  of  Puget  Sound,  model 
predictions  are  much  less  accurate.  Similar  behavior  is  seen  in  Prince  William 
Sound,  AK.  These  areas  are  ideal  candidates  for  regional-scale  high-resolution 
models  which  could  provide  a  more  precise  representation  of  the  bathymetry 
and  coastal  geometry  and  be  forced  with  the  results  from  the  coupled  global/ 
ENPACT  model. 

The  results  obtained  from  this  study  are  generally  less  accurate  than  those 
obtained  from  a  similar  modeling  effort  along  the  U.S.  east  coast,  the  Gulf  of 
Mexico,  and  the  Caribbean  Sea  (Westerink,  Luettich,  and  Scheffher  1993). 
While  it  is  not  entirely  clear  why  this  should  be,  the  most  likely  explanations 
would  seem  to  be  inaccurate  bathymetry  and  poor  offshore  boundary  condi¬ 
tions.  It  is  also  possible  that  the  model  needs  to  be  run  with  a  radiative  boun¬ 
dary  condition  to  allow  energy  to  propagate  out  of  the  domain.  This  may  be 
particularly  important  for  the  eastern  North  Pacific,  which  has  much  steeper 
(and  hence  more  reflective)  bottom  topography  and  much  narrower  (and  hence 
less  dissipative)  shelves  than  along  the  U.S.  east  coast. 
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The  ENPACT  model  is  a  large-domain  coastal  ocean  model  with  a  geo¬ 
metrically  and  hydrodynamically  simple  open  boundary.  The  large  domain 
simplifies  open  boundary  specification  and  allows  inclusion  of  large-scale 
effects  arising  from  the  Gulf  of  Califomia  and  the  Guatemala  and  Panama 
Basins.  The  inclusion  of  these  large-scale  effects  adds  little  to  the  size  of  the 
discrete  problem  since  the  model  uses  a  variable-spaced  grid. 

The  model  was  forced  with  five  astronomical  constituents  on  the  open 
boundary  derived  from  a  state-of-the-art  global  tidal  model.  In  addition,  a 
realistic  tidal  potential  forcing  was  applied  to  the  interior  of  the  domain.  The 
simulations  were  intended  to  be  entirely  predictive,  so  no  tuning  or  calibration 
procedures  were  performed. 

Comparison  between  model  results  and  field  data  was  made  at  72  stations 
distributed  throughout  the  domain  for  the  Kj,  Oj,  Mj,  Sj,  and  Nj  constituents. 
Average  amplitude  errors  for  aU  the  stations  range  from  15.3  to  65.5  percent, 
and  average  errors  in  phase  range  ft-om  5.0  to  66.2  deg.  Some  modded  con¬ 
stituents  (e.g.,  Oj)  show  excellent  agreement  with  observations,  while  others 
(e.g.,  Sj)  are  poorly  predicted.  Phase  errors  along  the  coast  are  quite  sensitive 
to  the  position  of  offshore  amphidromes.  Even  in  the  deep  ocean,  there  is  no 
close  match  between  predictions  and  observations.  Potential  reasons  for  these 
discrepancies  include  inaccurate  bathymetry,  poor  elevation  open  boundary 
conditions,  and  the  lack  of  a  radiative  condition  along  the  model  open 
boundary. 
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Alaska 


Figure  1.  The  Eastern  North  Pacific  Tidal  (ENPACT)  model  computational  domain;  approximate  locations  of  known 
amphidromes  determined  from  Qeosat  altimetry  (Cartwright,  Ray,  and  Sanchez  1991)  are  denoted  by  a 
triangle  and  the  constituent 


Figure  2.  Bathymetry  contours  (meters)  within  the  ENPACT  model  domain 


Graded  finite  element  grid  with  resolution  ranging  from  0.5  deg  to  less  than  7.5  min 
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Figure  6.  Computed  contours  for  elevation  amplitudes  (meters)  and  phases  (degrees  relative  to  Greenwich  Mean  Time) 
for  the  two  diurnal  astronomical  tidal  constituents  used  for  boundary  and  interior  domain  forcing. 
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Figure  7.  Computed  contours  for  elevation  amplitudes  (meters)  and  phases  (degrees  relative  to  Greenwich  Mean  Time) 
for  the  three  semidiurnal  astronomical  tidal  constituents  used  for  boundary  and  interior  domain  forcing. 
(Sheet  1  of  6) 
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Figure  8.  Time  histories  of  surface  elevations  at  representative  stations  within  the  ENPACT 
model  domain.  These  plots  have  been  constructed  using  the  Ki,  Oi,  M2,  S2,  and  N2 
tidal  constants  from  the  ENPACT  model  (solid  line)  and  from  observed  data  (dashed 
line).  (Sheet  1  of  10) 


ELEVATION  (meters) 


ELEVATION  (meters) 


ELEVATION  (meters) 


ELEVATION  (meters) 


Figure  8.  (Sheet  5  of  10) 
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Table  1 

Observational  Tidal  Elevation  Station  Locations 
and  Record  Durations 

station 

Number 

Station  Name 

Position 

(longitude  -  latitude) 

Location 

Record 

Length 

1 

Callao,  Peru 

coast 

369  days 

2 

Port  Liberty,  Equador 

coast 

358  days 

3 

Buenaventura,  Columbia 

coast 

369  days 

4 

Puerto  Armuelles,  Panama 

coast 

369  days 

5 

Puntarenas,  Costa  Rica 

coast 

369  days 

6 

San  Juan  Del  Sur,  Nicaragua 

coast 

369  days 

7 

La  Union  (Cutuco),  El  Salvador 

coast 

369  days 

8 

San  Jose.  Guatemala 

coast 

369  days 

9 

Salina  Cruz,  Mexico 

coast 

369  days 

10 

Acapulco,  Mexico 

coast 

369  days 

11 

Isla  Socorro,  Mexico 

coast 

369  days 

12 

Manzanillo,  Mexico 

coast 

369  days 

13 

Puerto  Vallarta,  Mexico 

coast 

369  days 

14 

Cabo  San  Lucas,  Mexico 

coast 

369  days 

15 

Mazatlan,  Mexico 

coast 

369  days 

16 

Topolobampo,  Mexico 

369  days 

17 

Guaymas,  Mexico 

coast 

369  days 

18  I 

Isla  Guadalupe,  Mexico 

coast 

369  days 

19 

Ensenada,  Mexico 

coast 

369  days 

20 

San  Diego,  CA 

coast 

2  years 

21 

Scripps  Institute,  La  Jolla,  CA 

coast 

369  days 

22 

Wilson  Cove,  San  Clemente  Is.,  CA 

coast 

29  days 

23 

Avalon,  Catalina  Island,  CA 

coast 

29  days 

24 

Ocean  Pier,  Balboa,  CA 

coast 

369  days 

25 

Outer  Harbor,  Los  Angeles,  CA 

coast 

3  years 

26 

Santa  Monica,  CA 

coast 

369  days 

27 

Port  Hueneme,  CA 

I^SI^QIIII^^QIIII 

coast 

369  days 

(continued)  | 

Source  of  Tidal  Data:  International  Hydrographic  Organization  Tidal  Constituent  Bank  (1990)  | 

Table  1  (Continued) 

station 

Number 

station  Name 

Position 

(iongitude  -  iatitude) 

Location 

Record 

Length 

28 

Avila,  CA 

coast 

369  days 

29 

Morro  Bay,  CA 

coast 

163  days 

30 

Monterey,  CA 

coast 

366  days 

31 

Golden  Gate,  San  Francisco,  CA 

coast 

11  years 

32 

lAPSO  30-2.1.14 

shelf 

58  days 

33 

Humboldt  Bay,  CA 

coast 

369  days 

34 

Crescent  City,  CA 

coast 

369  days 

35 

Brookings,  OR 

coast 

369  days 

36 

Marshfield,  Coos  Bay,  OR 

coast 

369  days 

37 

South  Beach,  Newport,  OR 

coast 

365  days 

38 

Youngs  Bay,  Astoria,  OR 

coast 

369  days 

39 

Toke  Point,  Willapa  Bay,  WA 

coast 

369  days 

40 

Neah  Bay,  WA 

coast 

369  days 

41 

Friday  Harbor,  WA 

coast 

365  days 

42 

Port  Renfrew,  Canada 

coast 

62  days 

43 

Cape  Scott,  Canada 

coast 

86  days 

44 

lAPSO  30-2.1.2 

shelf 

18  days 

45 

Wadhams,  Canada 

coast 

69  days 

46 

Milne  Island,  Canada 

coast 

50  days 

47 

Port  Simpson,  Canada 

coast 

69  days 

48 

Queen  Charlotte,  Canada 

coast 

62  days 

49 

Tasu  Sound,  Canada 

coast 

86  days 

50 

QIawdzeet  Anchorage,  Canada 

shelf 

87  days 

51 

Langara  Island,  Canada 

coast  ! 

39  days 

52 

Ketchikan,  AK 

coast 

3  years 

53 

Sitka,  AK 

coast 

2  years 

54 

Yukutat,  AK 

coast 

369  days 

55 

Cape  Whitshed,  Orca  Inlet,  AK 

coast 

29  days 

56 

Orca  Bay,  Prince  William  Sound,  AK 

coast 

29  days 

57 

Middleton  Island,  AK 

coast 

29  days 

(continued) 

Table  1  (Concluded) 

Station 

Number 

Station  Name 

Position 

(longitude  -  latitude) 

Location 

Record 

Length 

shelf 

145  days 

-  58 

I APSO  30-2.1.6 

59 

lAPSO  30-2.1.5 

shelf 

170  days 

60 

Seldovia,  AK 

coast 

365  days 

61 

Shaw  Island,  Cook  Inlet,  AK 

coast 

29  days 

62 

Kodiak,  AK 

coast 

369  days 

63 

King  Cove,  AK 

coast 

239  days 

64 

Peterson  Bay,  Sanak  Island,  AK 

162_38'  W  -  54_23'  N 

coast 

15  days 

65 

Battra  Island,  Galapagos  islands 

coast 

366  days 

66 

1  APSO  30-2.1.13 

124_.26'  W  -  27_45‘  N 

deep 

ocean 

26  days 

67 

lAPSO  30-2.1.11 

119_48'W  -  31_2‘N 

deep 

ocean 

37  days 

68 

lAPSO  30-2.1.12 

120_51‘W  -  32_14'N 

deep 

ocean 

8  days 

69 

lAPSO  30-2.1.9 

145^0’  W  -  34_0‘  N 

deep 

ocean 

15  days 

70 

lAPSO  30-2.1.1  (Cobb  Seamount) 

130^49'  W  -  46_,46'  N 

deep 

ocean 

203  days 

71 

Union  Seamount 

132_47‘  W  -  49_^35'  N 

deep 

ocean 

10  days 

72 

lAPSO  30-2.1.3 

145_43‘W  -  58__46'N 

deep 

ocean 

40  days 

Table  2 

Range  of  Amplitude  Values  for  Open  Boundary  Forcing 
Obtained  from  a  Globai  Ocean  Model 
(Le  Provost,  et  al.,  1994) 


Species,  / 

n 

Constituent 

Maximum  Ampjn ,  m 

Minimum  Ampjn ,  m 

1 

1 

K«|,  luni-solar 

0.4425 

0.0194 

2 

Oi,  principal  lunar 

0.2974 

0.0063 

2 

1 

M2,  principal  lunar 

0.6021 

0.0932 

2 

S2,  principal  solar 

0.1671 

0.0319 

3 

N2,  elliptical  lunar 

0.1218 

0.0226 

Table  3 

Tidal  Potential  Constants  for  the  Principal  Tidal  Constituents 
(from  Reid,  1990)  and  the  Associated  Effective  Earth  Elasticity  Factor 
(from  Wahr,  1981) 


Species,  / 

n 

Constituent 

Tjnt  hr 

Cjn.m 

o-jn 

1 

1 

Ki,  luni-solar 

23.934470 

0.141565 

0.736 

2 

Oi,  principal  lunar 

^.819342 

0.100514 

0.695 

2 

1 

M2,  principal  lunar 

12.420601 

0.242334 

0.693 

2 

S2,  principal  solar 

12.000000 

0.112841 

0.693 

3 

N2,  elliptical  lunar 

12.658348 

0.046398 

0.693 

Table  4 

Comparison  of  Amplitude  Errors  Between  Calcuiated  and  Observed 
Values,  for  Constituent  /,  in  Various  Regions 

of  the  Domain 


C  calc  -oba 


(see  Equation  11) 


Constituent 

Entire 

Domain 

Stations 

1-72 

South  and 
Central 
America 
Stations 

1-19 

U.S.  West 
Coast 
Stations 
20-41 

Canada  and 
Alaska 
Stations 
42-64 

Offshore 

Stations 

65-72 

Kl 

0.153 

0.064 

0.237 

0.085 

0.028 

Ol 

0.165 

0.088 

0.228 

0.133 

0.033 

M2 

0.227 

0.362 

0.338 

0.145 

0.215 

S2 

0.655 

0.535 

1.225 

0.531 

0.687 

N2 

0.184 

0,299 

0.251 

0.118 

0.149 

Table  5 

Comparison  of  Phase  Errors  Between  Calculated  and 
Observed  Values,  E“phatf  •  for  Constituent  y,  in  Various  Regions 

of  the  Domain 


•  caic-obs 
phase 


(see  Equation  1 2) 


Constituent 


K1 


01 


M2 


S2 


N2 


Entire 

Domain 

Stations 

1-72 


South  and 

Central 

America 

Stations 

1-19 


U.S.  West 
Coast 
Stations 
20-41 


Canada  and 
Alaska 
Stations 
42-64 

Offshore 

Stations 

65-72 

6.5 

3.7 

6.8 

8.0 

26.4 

38.7 

40.6 

94.6 

30.0 

41.4 
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Appendix  A 

Resynthesis  of  Time  Histories 
Using  the  Cotidal  Charts 


Resynthesized  time  histories  are  useful  for  examining  the  general  structure 
of  tides  and  for  overall  comparisons  between  model  results  and  observed  data. 
In  order  to  resynthesize  a  time  history  from  the  cotidal  charts  presented  in 
Figures  6  and  7,  the  following  steps  can  be  applied; 


Step  1 

Interpolate  amplitudes  and  phases  for  the  five  astronomical  constituents 
presented  in  this  report  using  fte  most  detailed  set  of  cotidal  charts  available. 
For  example,  in  order  to  obtain  values  for  Cape  Whitshed,  Orca  Inlet,  AK 
(located  at  145°  55’  W  and  60°  28’  N),  the  set  of  cotidal  charts  presented  in 
Figures  6  and  7  are  interpolated.  The  adjacent  amplitudes  A,  and  'F  phases 
can  be  readily  obtained; 


Constituent 

Aj,  in 

S',,  deg 

K, 

0.504 

273 

o, 

0.323 

256 

M2 

1.360 

307 

Sg 

0.619 

352 

N2 

0.252 

287 

step  2 

Surface  elevations  as  a  fimction  of  time  can  now  be  computed  using  the 
following  relationship; 
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A1 


(Al) 


cos[  m.  (t  -  +  V .  (tj  -  (X,  3')] 

where 

C(^.  >’•  0  =  surface  elevation  at  the  site  of  interest  as  a  function  of 
time 

Ai(x,  y)  =  amplitude  for  the  t'*  constituent 

*p.  (x,  y)  =  phase  for  the  i'*  constituent  (in  radians  per  second) 

f^tg)  =  nodal  factor  relative  to  reference  time  tg  for  the 
i**  constituent 

v^tg)  =  equilibrium  argument  relative  to  reference  time  tg  for  the 
constituent 


For  the  five  constituents  of  interest,  cO;  has  the  following  values: 


Constituent 

(radians/sec) 

K, 

0.000072921 158358 

Oi 

0.000067597744151 

Ma 

0.000140518902509 

Sa 

0.000145444104333 

Nj, 

0.000137879699487 

Guidelines  as  to  how  to  compute /f^and  vf^are  given  by  Schureman  (1941) 
and  Foreman  (1977).  If  equilibrium  tides  are  computed, /f^and  vf  ^should  be 
specified  as  1.0  and  0.0  for  all  constituents.  A  resynthesized  signal  for  the 
data  is  presented  between  30  and  40  days  in  Figure  Al. 


Step  2  (Alternate) 

As  an  alternative  to  computing  a  time  series  according  to  Equation  Al,  the 
constiment  tide  record  generation  program  in  the  Automated  Coastal  Engineer¬ 
ing  System  (ACES)  (Leenknecht,  Szuwalski,  and  Sherlock  1990)  can  be  used. 
This  program  requires  constituent  data  for  the  specified  location  to  predict  tidal 
elevations  at  that  locale  for  a  given  time  period.  The  tidal  constituents 
contained  in  this  report  can  be  used  as  input  to  the  ACES  program  with  one 
modification.  The  Greenwich  epochs  4*  shown  in  Equation  Al  must 


A2 
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55  Cape  Whitshed,  Orca  Inlet,  AK 


Figure  A1 .  Example  time  history  of  resynthesized  surface  elevations  from  the 
K,.  O,,  Mg.  Sg,  and  Ng  tidal  constituents  at  a  station  in  the  Eastern 
North  Pacific. 


be  converted  to  local  epochs  k.  This  transformation  is  according  to  the  fol¬ 
lowing  relationship: 

Greenwich  epoch  \|r-  =  k-  +  pL  (A2) 


where  k  is  the  local  epoch  needed  for  input  to  the  ACES  program  to  recon¬ 
struct  a  local  tide,  p  is  the  coefficient  of  the  constituent  (i.e.,  2  for  Mj,  1  for 
K^,  etc.),  and  L  is  the  longitude  of  the  location  for  which  a  tidal  prediction  is 
desired.  By  making  the  above  adjustment,  this  report  can  be  used  in  conjunc¬ 
tion  with  the  ACES  program  to  yield  tidal  predictions  at  any  location  along  the 
west  coast  of  the  United  States,  Canada,  and  Central  and  part  of  South 
America. 


Step  3  (Optional) 

Comparisons  with  observational  data  are  often  useful  in  evaluating  model 
accuracy.  Observational  data  available  from  the  International  Hydrographic 
Organization  Tidal  Constituent  Bank  and  other  organizations  are  sometimes 
referenced  to  the  epoch  of  the  local  station.  The  following  transformation 
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A3 


converts  these  data  to  the  Greenwich  epoch  to  allow  comparisons  with 
ADCIRC  results: 


Greenwich  epoch  \|r.  =  g-  +  a- 


15 


(A3) 


where  gi  is  the  station’s  local  epoch,  a-  is  the  speed  of  the  constituent  argu¬ 
ment  in  degrees/hour,  and  S  is  longitude  of  the  station’s  time  zone  meridian, 
positive  for  west  and  negative  for  east  (Schureman  1941). 


Constituent 

Sf  (deg/hour) 

K, 

15.0410686 

o, 

13.9430356 

M2 

28.9841042 

Sa 

30.0000000 

N2 

28.4397295 

A4 
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